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Abstract. Studies of stellar magnetism at the pre-main sequence phase can provide important 
new insights into the detailed physics of the late stages of star formation, and into the observed 
properties of main sequence stars. This is especially true at intermediate stellar masses, where 
magnetic fields are strong and globally organised, and therefore most amenable to direct study. 
This talk reviews recent high-precision ESPaDOnS observations of pre-main sequence Herbig 
Ae-Be stars, which are yielding qualitatively new information about intermediate-mass stars: 
the origin and evolution of their magnetic fields, the role of magnetic fields in generating their 
spectroscopic activity and in mediating accretion in their late formative stages, and the factors 
influencing their rotational angular momentum. 
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1. Introduction 

1.1. Magnetism and rotation in the main sequence A and B stars 

Between about 1.5 and 10 M0, at spectral types A and B, about 5 % of main sequence 
(MS) stars have magnetic fields with characteristic strengths of about IkG. Such stars 
also show important chemical peculiarities and are thus usually called the magnetic chem- 
ically peculiar Ap/Bp stars. The strength of the magnetic fields of these stars cannot be 
explained by an envelope dynamo as in the sun. Until now, the most reliable hypothe- 
sis has been to assume a fossil origin for these magnetic fields. This hypothesis implies 
that the stellar magnetic fields are relics from the field present in the parental interstellar 
cloud. Its also implies that magnetic fields can (at least partially) survive the violent phe- 
nomena accompanying the birth of stars, and can also remain throughout their evolution 
and until at least the end of the MS, without regeneration. 

According to the fossil field model, we should observe magnetic fields in some pre-main 
sequence (PMS) stars of intermediate mass, the so-called Herbig Ae/Be stars. How- 
ever no magnetic fi eld was observed up to recently in these stars (except HD 104237, 
Donati et al. |[l997h . Can we obtain some observational evidence of the presence of mag- 
netic fields during the PMS phase of evolution, as predicted by the fossil field hypothesis? 
If some Herbig Ae/Be stars are discovered to have magnetic fields, is the fraction of mag- 
netic to non-magnetic Herbig stars the same as the fraction for main sequence stars? Is 
the magnetic field in Herbig stars strong enough to explain the strength of that of Ap/Bp 
stars? 
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Chemical peculiarities and magnetism are not the only characteristic properties ob- 
served in the Ap/Bp stars. Most magnetic MS stars have rotation periods (typically of 
a few days) that are several times longer than the rotation periods of non-magnetic MS 
stars (a few hours to one day). It is usually believed that magnetic braking, in particular 
during PMS evolution, when the star can exchang e angular momentum with its massive 
accre tion disk, is responsible for this low rotation ( Stepien 2000t St^pieh fc Landstreet I 



20021) . An alternative involves a rapid dissipation of the magnetic field during the early 



stages of PMS evolution for the fastest rotators, due to strong turbulence induced by 
rotational shear develo ped under the surface of the stars, as the convection do in the 
solar-type stars (see e.g. Lignieres et al. 19961 ). In this scenario, only slow rotators could 
retain their initial magnetic fields, and evolve as magnetic stars to the main sequence. So 
the question to be addressed is the following: does the magnetic field control the rotation 
of the star, or else does the rotation of the star control the magnetic field? We propose 
that this question can be answered by studying rotation and magnetic fields in Herbig 
Ae/Be stars. 



1.2. The Herbig Ae/Be stars 



The Herbig Ae/Be stars are intermediate-mass pre-main sequence stars, and therefore 
the evolutionary progenitors of the MS A and B stars. They are distinguished from the 
classical Be stars by their IR emission and the association with nebulae, characteristics 
which are due to their young age. 

They display many observational phenomena often associated wit h magnetic activity . 
First, high ionised lin es are observed in the spectra of some stars (e.g. Bouret et al. Ill997 : 
iRoberge e t al. '^'2001'), and X-ray emission have been detected, coming from some Herbig 
stars (e.g. Hamaguchi et al. 2005). In active cool stars, many of these phenomena are 
produced in hot chromospheres or coronae. Some authors mentioned rotational modula- 
tion of resonance lines which they s peculate m ay be due to rotation modulation of winds 
structured by magnetic field (Prade rie et al. 1986; Catala et al. 1989, 1991, 1989). 

In the literature we find many clues of the presence of circumstellar disks around 
these stars, f rom spectroscopic data showing s trong emission, and also from photomet- 
ric data (e.g. Mannings fc Sargent 1997 . 2000l ). Recently, using coronagraphic data and 
interferometric dat a, some authors have also found direct evide nce of circumstellar disks 
around these stars (jGradv et al. 1 ll999ll2000l :l lEisner et al. 1 l2003h . A careful study of these 
disks shows that th ey have similar properties to the disk of their low mass counterpart 
( Natta et al. I I2OOII) . the T T a,uri stars, whose the emission lines are exp l ained by mag- 
netospheric accre tion models ( Koenigl et al. 1991; MuzeroUe et al. 1998 . 200 ll ). Finally 



Muzerolle et al. 



([2004) have sucessfully applied their magnetospheric accretion model 
to Herbig stars to explain the emission lines in their spectra. 

For all these reasons we suspect that the Herbig stars may host large-scale magnetic 
fields that should be detectable with current in strumentation. However, many authors 



1997t iHubrig et al. 



tried to detect such fields without muc h success (ICatala et al 
200I IWade et al" 



19931 . ll989l : [Donati et al 



20071 ). But in 2005, a new high-resolution spec- 



tropolarimeter, ESPaDOnS, has been installed at the canada-france-hawaii telescope 
(CFHT). We therefore decided to proceed to survey many Herbig stars in order to inves- 
tigate rotation and magnetism in the pre-main sequence stars of intermediate mass. 
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Figure 1. Herbig stars plotted in a HR diagram. The red squares are the magnetic Herbig stars. 
The PMS evolutionary tracks are plotted for different masses (full lines). The birthline for 10~^ 
and 10~* MQ.yr"^ mass accretion rate are plotted in dashed line l|Palla fc Stahler1ll993l ). The 
dash-dotted line is the ZAMS. 

2. Observations and reduction 

2.1. Our sample 



We ha ve selected Herbig stars in the catalogues from I The et al. I (|1994D and lVieira et al. 



with a visual magnitude brighter than 12. Our sample contain 55 stars which 
have masses ranging from 1.5 M0 to 15 Mq with all ages between the birthline and the 
zero-age main sequence (ZAMS). In Fig. [1] are plotted the stars of our sample in an HR 



diagr am, as well as the evolutionary t racks computed with t he CESAM code ([Morel 



1997 ). and the birthlines computed bv IPalla &: Stahler" (1993) with two mass accretion 
rates during the protostellar phase : 10~^ M0.yr~^ and 10"'' M0.yr~^. 



2.2. 



Observations and reduction 

Our data were obtained using the high resolution spectropolarimeter ESPaDOnS installed 
on the 3.6 m Canada- France-Hawaii Telescope (Donati et al. 2007, in preparation) during 
many scientific runs. 

We used this instrument in polarimetric mode, generating spectra of 65000 resolution. 
Each exposure was divided in 4 sub-exposur es of equal time in order to compute the 
optimal extraction of the polarisation spectra ([Donati et al. I I1997I , Donati et al. 2007, in 
prep.). We recorded only circular polarisation, as the Zeeman signature expected in linear 
polarisation is about one order of magnitude lower than circular polarisation. The data 
were reduced using the "Libre ESpRIT" package especially developed for ESPaDOnS, 
and installed at the CFHT (Donati et al. 1997, Donati et al. 2007, in preparation). After 
reduction, we obtained the intensity Stokes / and the circular polarisation Stokes V 
spectra of the stars observed. 



We then applied the Least Squares Deconvolution procedure to all spectra ([Donati et al 



19971 ). in order to increase our signal to noise ratio. This method assumes that all lines of 
the intensity spectrum have a profile of similar shape. Hence, this supposes that all lines 
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Figure 2. Stokes / (bottom) and Stokes V (up) LSD profiles plotted for the 4 magnetic stars 
(right and middle) and two undetected stars (left). Note the amplification factor in V. 



are broadened in the same way. We can therefore consider that the observed spectrum is 
a convolution between a profile (which is the same for all lines) and a mask including all 
lines of the spectrum. We therefore apply a deconvolution to the observed spectrum using 
the pre-computed mask, in order to obtain the average photospheric profiles of Stokes I 
and V. In this procedure, each line is weighted by its signal to noise ratio, its depth in the 
unbroadened model and its Lande factor. For each star we used a mask computed using 
"extract stellar" line lists obtained from the Vienna Atomic Line Database (VALDjj, 
with effective temperatures and log 5 suitable for each star (Wade, Alecian et al. 2007, in 
prep.). We excluded from this mask hydrogen Balmer lines, strong resonance lines, lines 
whose Lande factor is unknown and emission lines. The results of this procedure are the 
mean Stokes / and Stokes V LSD profiles (Fig. [2])- 



3. Results 

3.1. Discovery oj magnetic fields in Herbig stars 
Thanks to the high performance of the instrument ESPaDOnS and t o the LSD method . 



we have discovered four new magn etic Herbig Ae/Be stars (|Wade et al. 112005 



Catala et alTlbOOeHAiecian et al. II2007I ). In Fig. [2] is plotted the Stokes / and V profiles 



of each new discovered magnetic Herbig Ae/Be star : HD 200775, HD 72106, V380 Ori 
and HD 190073, and of four stars in which a magnetic field has not been detected (the 
undetected stars, hereafter). Contrary to the undetected stars, the Stokes V profiles of 
the magnetic stars are not null and display a strong Zeeman signature, of the same width 
of the photospheric / profile, characteristic of the presence of a magnetic field in the stars. 

These 4 detections among our sample of 55 stars lead to the conclusion that 7% of 
Herbig Ae/Be stars are magnetic. The projection of the distribution of magnetic and 
non-magnetic main sequence A and B stars on the pre-main sequence phase, assuming 
a fossil field hypothesis, predict that between 5 and 10 % of Herbig stars should be 



f http://www.astro.univie.ac.at/~vald/ 
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Table 1. Fundamental, geometrical and magnetic parameters of the magnetic Herbig Ae/Be 

stars. References : 1 : [Alecian et a l. (2007), 2 : Folsom et al. (2007), in prep., 3 : Alecian et al. 

(2007), in prep., 4 : ICatala et aTl lj200ll ). 



Star 


S.T. 


V sin i 


age 


P 


Bp 




i n 


ddip 


-Bp(ZAMS) 


Ref. 






(km.s^^) 


(Myr) 


(d) 


(kG) 






R* 


(kG) 




HD 200775 


B2 


26 


0.1 


4.328 


1 


78 


13 


0.1 


3.6 


1 


HD 72106 


AO 


41 


10 


0.63995 


1.5 


58 


23 





1.5 


2 


V380 oriX 


A2 


9.8 


2.8 


[7.6,9.8] 


1.4 


[90, 85] 


[36, 49] 





2.4 


3 


HD 190073^ 


A2 


8.5 


1.5 




[0.1,1] 


[0,90] 


[0,90] 




[0.3,3] 


4 



f Work in progress. We need more data of V380 Ori to choose between the 7.6 and the 9.8 
periods. Therefore two solutions are possible for /3 and i 
J Although we have observed HD 90073 over more that 2 years, no variation of the Stokes V 

profile has yet been detected. 



magnetic, which is consistent with our observations. We therefore bring a new strong 
argument in favoured of the fossil field hypothesis (Wade, Alecian et al. 2007, in prep.). 

3.2. Topology and intensity of the magnetic fields 

3.2.1. The magnetic Herbig stars 

We determined the topology and the intensity of the magnetic fields of the four mag- 
netic Herbig stars in order to compare them to the magne tic M S A an d B stars. With 



this aim, we used the oblique rotator model described by IStift I (|l975l ) . We consider a 
dipole placed at a distance ddip on the magnetic axis of a spherical rotating star with a 
magnetic intensity at the pole Bp. The rotation axis of the star is inclined at an angle i 
with respect to the line of sight and makes an angle f3 with the magnetic axis. 



According to lLandi deglTnnocenti I (|l973l ). in the weak field approximation, the Stokes 



V profile is proportional to the magnetic field projected onto the line of sight and in- 
tegrated over the surface of the star {Bg, the longitudinal magnetic field, hereafter). As 
the star rotates, the visible magnetic changes, resulting in variation of Bi. Therefore the 
Stokes V profile changes with the rotation phase. 

In order to determine the geometrical and magnetic parameters i, [3, Bp and ddip, 
as well as the rotation period P of the star, we observed the stars at different rotation 
phases and fit simultaneously all the Stokes V profiles observed for each star. With 
this aim we calculated a grid of V profiles, using the oblique rotator model, for each 
date of observations, varying the five parameters. Then, for each star, we applied a 
minimisation to find the best model which matches simultaneously all the V profiles 
observed. Fig. [3] shows the result of our fitting procedure fo r one star : HD 2007 75. The 



synthetic Stokes V are superimposed on the observed ones ( Alecian et al. 2007t ). 

The values of the geometrical and magnetic parameters are summarized in Table [1] for 
each stars. In the case of HD 190073, the topology and the intensity of its magnetic field 
are not constrained, because, during 2 years observations, the Stokes V profile has not 
been observed to vary. There are 3 possible explanations for this : the inclination i is very 
small, the obliquity angle /3 is very small, or the rotation period of the star is very long. 
More observations will allow us to discard two of these solutions. However the stability 
of the magnetic field over more than 2 years and the shape of the Stokes V profiles lead 



us to the conclusion that this star hosts a large-scale fossil magnetic field (jCatala et al. 



The success of our fitting procedure for the 3 stars HD 200775, HD 72106 and V380 
ori, as well as our discussion on HD 190073, lead to the conclusion that the magnetic 
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Figure 3. Stokes V profiles of HD 200775 superimposed to the syntlietic ones, corresponding 
to our best fit. The rotation p hase and the error b ars are indicated on the left and on the right 
of eache profiles, respectively (jAlecian et al. II2007I ). 

Herbig Ae/Be stars host globally dipolar magnetic fields, similar to the Ap/Bp 
stars. 

Assuming the conservation of magnetic flux during the PMS evolution, and using the 
radius of the stars and their predicted radius on the ZAMS for the same mass, we can 
estimate the magnetic intensity at their surface they will have when they will reach 
the ZAMS (see Tabled]). We found intensities ranging from 300 G to 3.6 kG, which is 
very close to what is observed in the Ap/Bp stars. Hence we bring 2 more new strong 
arguments in favour of the fossil fleld hypothesis. 

3.2.2. The undetected stars 

In the case of the other 51 stars no magnetic field was detected. However, a non- 
detection does not mean that the star does not host a magnetic field. The star may host 
a magnetic field which is too weak to be detected with our method, or we way observed it 
at a rotation phase where the Stokes V profile was too weak to be detected. However, for 
two of these stars (HD 139614, HD 169142), with wsini similar to the magnetic Herbig 
stars, we obtained data during many successive nights, with a 3ct uncertainty on the 
measured Bi less than 50 G, which is typical of the uncertainties on Be obtained in the 
magnetic Herbig stars. Therefore, if these stars would host a magnetic field similar to 
the magnetic Herbig stars, we would have detect it. This implies that if these stars host 
a magnetic field, their magnetic intensities are much lower than those of the magnetic 
Herbig stars. We will consider only these two well-constrained stars in the following, 
because more work need to be done to give conclusions on the magnetic intensity of the 
other stars. 



4. Consequences on magnetospheric accretion in Herbig stars 

We can use these intensities to test if magnetospheri c accretion can occur in these stars. 



Following IWade et al. (120071) we used th e models of iKoenigl et al. I ( 199ll ). [Shu et al. 



(1994) and lCollier Cameron fc Campbell I (11993 ). which consider a dipole magnetic field 
aligned with the rotation axis, and coupled with the star until the corotation radius or a 



Magnetism, rotation and accretion in Herbig Ae-Be stars 



125 



fraction of the corotation radius. The matter from the disk is therefore funnelled along 
the magnetic lines until the surface of the star. We used the equations, given by these 
authors, of the polar magnetic field in function of the parameters of the stars to calculate 
the minimum intensity requ ired to trigger magnetospheric accretion in these stars (see 



also Jons-Krull et al. 1999f) . We considered a mass accretion rat e of 10 * M, 



which is the typical value in Herbig stars (e.g. Blondel et al. 19931 ) 



0-yr" 



In the case of t he magnetic stars we found for the models of iKoenigl et al. I (|1991I ) and 



Shu et al. I (|1994| ). that only one stars on four have a stron g enough magnetic intensity to 



cause magnetospheric accretion, whereas with the model of lCollier Cameron &: Campbell 
( 1993h . the four stars can have magnetospheric accretion. 

In the case of the two well constrained undetected stars, considering that they could 
host a magnetic field with intensity lower than few 100 G, the minimum field intensity 
required is too high to have magnetospheric accretion. 

Therefore, according to the models that we considered, and taking into account our 
observations, magnetospheric accretion cannot occur in all Herbig stars. 



5. Conclusions 

We used the new spectropolarimeter ESPaDOnS installed at the CFHT to proceed in 
a survey of the Herbig stars, in order to investigate their rotation and magnetic field. 
We discovered four magnetic stars whose field topology is similar to the MS magnetic 
A-B stars. We also show that the magnetic intensities of these fields and the proportion 
of magnetic Herbig stars can explain the magnetic intensity and the proportion of mag- 
netic fields among the MS stars, in the context of fossil field model. We therefore bring 
fundamental arguments in favour of this hypothesis. 

The four magnetic Herbig stars are slow rotators (wsini < 41 km.s^^) which supports 
that magnetic Herbig Ae/Be stars are the progenitors of the magnetic Ap/Bp stars. 
Among these magnetic stars two have very low vsini (< 10 km.s~^) and are very young 
(age< 2.8 Myr). Assuming these stars are true slow rotators, this implies that there exists 
a braking mechanism which acts very early during the PMS evolution of the intermediate 
mass stars. We could think that this braking mechanism has a magnetic origin, although 
among the undetected stars we also observe same stars with small vs'mi (^ 15 km.s^^). 
The nature of the braking mechanism requires addition study. 

Finally it has been proposed that all Herbig stars should undergo magnetospheric ac- 
cretion, as the spectra of all stars show similar emission. However, we calculated the min- 
imum polar magnetic intensity of the magnetic Herbig stars and of two well-constrained 
undetected stars to get magnetospher ic accretion, using three different models which have 



been developed for the T Tauri stars (jKoenigl et al. II1991I : ICoUier Cameron fc Campbell 



1993t IShu et al. Ill994f ). Considering the intensity of the magnetic stars, as well as the 
maximum magnetic intensity of the two well-constrained undetected stars, if they host 
a magnetic field, our first conclusion is that magnetospheric accretion cannot occur in 
all the Herbig stars (a statistic study taking into account all the undetected stars is in 
progress in order to confirm this result). Therefore, cither the models that we used are 
not well adapted to the Herbig stars, or the emission lines are not only produced by mag- 
netospheric accretion. A thorough observational and theoretical study of the emission in 
the spectra, as well as the surroundings of the Herbig Ae/Be stars is necessary to better 
understand the interaction of these stars with their surroundings. 
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